A df particle method for gyrokinetic simulations with kinetic electrons and electromagnetic perturbations

A df particle simulation method is developed for solving the gyrokinetic-Maxwell system of equations that describes turbulence and anomalous transport in toroidally confined plasmas. A generalized split-weight scheme is used to overcome the constraint on the time step due to fast parallel motion of the electrons. The inaccuracy problem at high plasma b is solved by using the same marker particle distribution as is used for df to evaluate the bmi=meAk term in Ampere s equation, which is solved iteratively. The algorithm is implemented in three-dimensional toroidal geometry using field-line-following coordinates. Also discussed is the implementation of electron–ion collisional effects which are important when kinetic electron physics is included. Linear benchmarks in toroidal geometry are presented for moderate b, that is, b 1, but bmi=me 1. Nonlinear simulation results with moderate b are also presented. 2003 Elsevier Science B.V. All rights reserved. PACS: 52.65.+z; 52.35.Mw